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It is known that a single product shock measure in some ol one-dimensional driven-diflusive 
systems with nearest-neighbor interactions might evolve in time quite similar to a random walker 
moving on a one-dimensional lattice with reflecting boundaries. The non-equilibrium steady-state of 
the system in this case can be written in terms of a linear superposition of such uncorrelated shocks. 
Equivalently, one can write the steady-state of this system using a matrix-product approach with 
two-dimensional matrices. In this paper we introduce an equilibrium two-dimensional one-transit 
walk model and find its partition function using a transfer matrix method. We will show that there 
is a direct connection between the partition functions of these two systems. We will explicitly show 
that in the steady-state the transfer matrix of the one-transit walk model is related to the matrix 
representation of the algebra of the driven-diffusive model through a similarity transformation. The 
physical quantities are also related through the same transformation. 
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I. INTRODUCTION 



During the last couple of years, the time evolution of a 
single product shock measure in one-dimensional driven- 
diffusive systems H, 11, 01 been studied in a series of 
papers 0, H, B 0, H 3 ■ It has been observed that by ap- 
plying some constraints on the microscopic reaction rates 
of the system the shock front might have simple random 
walk dynamics while reflecting from the boundaries of the 
lattice. In this case the steady-state of the system can be 
easily obtained by considering a linear superposition of 
such product measures. Three families of two-state sys- 
tems with nearest-neighbor interactions have been intro- 
duced and studied in Q. This has been also generalized 
to the multi-species systems with long range interactions 
and even the systems with discrete time updating scheme 

The study of the steady-state properties of one- 
dimensional driven-diffusive systems has a long history. 
Different approaches have been used in order to find the 
stationary-state probability distribution function of these 
systems ranging from the Bethe ansatz to the matrix- 
product approach. In the latter scenario one writes the 
stationary-state of the system as a matrix element of 
product of non-commuting operators associated with dif- 
ferent states of each lattice site of the system. The alge- 
braic relation between these operators, sometimes called 
the algebra of the system, might have finite- or infinite- 
dimensional matrix representations pdj . In an attempt 
to understand the nature of the matrix representations of 
these algebras, it has been observed that for the systems 
in which a single product shock measure has a simple 
random walk dynamics under the time evolution gen- 
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erated by the so-called Hamiltonian of the system, the 
stationary-state probability distribution function can be 
expressed by a two-dimensional matrix representation. 
These matrices contain all of the information about the 
hopping rates of the shock front and also the densities of 
the particles on the left and the right hand side of the 
shock position. It has been confirmed that the conditions 
under which the product shock measure has a simple ran- 
dom walk dynamics are exactly those for the existence of 
a two-dimensional matrix representation for the algebra 
of the system [12|. 

Following the paper by Arndt [H investigations have 
shown that some of the concepts used in the equilib- 
rium statistical mechanics can be extended to the out-of- 
equilibrium systems. These concepts include the Yang- 
Lee description of phase transitions in equilibrium sys- 
tems. It is known that by defining an ad hoc partition 
function in terms of the microscopic reaction rates of 
an out-of-equilibrium system, one can apply the Yang- 
Lee approach to spot the transition points. On the 
other hand, in [3, [3, [iBl the connection between the 
steady-state of a one-dimensional driven-diffusive system 
i.e. the Asymmetric Simple Exclusion Process (ASEP) in 
random sequential updating scheme with that of a one- 
transit walk model was investigated. The investigations 
have been also extended to the case of ASEP with par- 
allel updating scheme in [Tt'I. 

The authors in ^15] have shown that by introducing two 
fugacities associated with the densities of contact points 
in the equilibrium one-transit walk model one can clearly 
explain the phases in the walk model. On the other hand, 
they have shown that if one replaces the two fugacities 
with the boundary rates of ASEP, a better understand- 
ing of the phase diagram of this non-equilibrium system 
can be obtained. In the present paper we extend the 
above mentioned ideas to the case where the shocks with 
simple random walk dynamics appear in the most gen- 
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eral one-dimensional driven-difFusive system. We define a 
one-dimensional driven-diffusive lattice model with open 
boundaries in which it is assumed that a single product 
shock measure can move on the lattice with random walk 
dynamics. One can introduce an ad hoc partition func- 
tion for this model in terms of its microscopic reaction 
rates. We also introduce an equilibrium one-transit walk 
model. As in pj| the partition function of this model 
can be calculated by defining proper fugacities. We aim 
to show that the physical quantities in both models are 
closely related. We believe that the connection between 
these two systems in the steady-state helps us understand 
the nature of applicability of equilibrium concepts to the 
non-equilibrium systems. 

This paper is organized as follows: We start with the defi- 
nition of the walk model and calculate some of its relevant 
physical quantities in terms of two proper fugacities. In 
order to find the partition function of this model we will 
use the transfer matrix method. Then we consider a gen- 
eral one-dimensional driven-diffusive system and show 
how its steady-state probability distribution function can 
be calculated using a matrix-product method and a two- 
dimensional matrix representation. We finally conclude 
by showing that the physical quantities of these two sys- 
tems are closely related if one relates the hopping rates 
of the shock front in the driven-diffusive system with the 
fugacities in the walk model. 

II. THE WALK MODEL 

Consider a random walker which moves on a two- 
dimensional path. Each path is defined on a diagonally- 
rotated square lattice. The paths start at (0, 0) and end 
at (2n,0). The walker can only move in the north-east 
(NE) or in the south-east (SE) direction. The walker can 
never take two consecutive steps to NE while it might 
take two consequent steps to SE only once. Each path 
either crosses the horizontal line only once (in this case 
the path can be factorized into two Dyck paths) or never 
crosses it (in this case the path is a single Dyck path). 
According to our definition there arc only ri + 1 different 
paths of this type. A typical path is illustrated in FIG. 
([1]). As can be seen, the height of each step can be 
or — 1. A fugacity zi is given to each down step and Z2 
to each up step, except those ending on the x-axis. The 
partition function of this model which contains 2n steps 
can be easily written as 

Zn{zi, Z2) = (ziZ2)"^n(2l, Z2) (1) 
I 

{s[,s'-,--- ,s'^^j,Si,Sj,--- ,s„} 



The matrix representations for the base vectors can be 
chosen as 

|l+> = |0+)= (J) , \l-) = \Q-)=[\y (6) 



in which 

''^ — n —n 

z„(zi,z2) -E-rv"^' = (2) 



The expression ([2]) has a simple interpretation given in 
[l6t . The weight associated with a path con- 

taining q contacts with the x-axis from above, is obtained 
by assigning a factor Zi^ to each contact with the x-axis 
from above, excluding an initial upward step, and a factor 
Z2^ to each contact from below, excluding a final upward 
step. As explained in T&l the prefactor (2:122)" does not 
change the critical behavior of the walk model. 
In what follows we will show that the partition function 
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FIG. 1: A typical one-transit walk. 



Zn{zi, 22) can be written as 

Z„(zi,Z2) = (i|r"ii?) (3) 

in which T is a two-step transfer matrix. We associate 
T° to each odd step and to each even step and write 
T — T°T'^ . In order to find the transfer matrix T and the 
vectors \R) and {L\ we first introduce four base vectors 
|si) = |1+), |S2) = |s'i) = |0+) and |4) = 10") with 
the following properties 

E|s,)(s,|=5]|s^)(s^|=X (4) 

in which I is an identity 2x2 matrix. Using the complete- 
ness of the base vectors the partition function Z„(2i, 22) 
can be expand as follows 



(5) 

I 

We assign these base vectors to the different vertices of 
each path according to the FIG. Now according to 



{L\s[) {s[\T°\s,){s,\T%s'Ms'^T°\s,){s, \ ■ ■ ■ 
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FIG. 2; Assigning the base vectors to the vertices 



the definition of the model only the following transitions 
are allowed 



(1+|T^|0+} = ^ 

(i-|r'=|o+) = o' 

(0+|T°|l+) = 1 
(0-|T°|l+} = 
(0+|i?,) = 1 
(L|0+) = 1 



(1+|T^|0-) = 
(l-|T«|0-) = ^ 

(o+|r°|i-) = 1 
(o-|r°|i-) = 1 

(L|o-) = o. 



(7) 



These equations determine the transfer matrices T° and 



rjiO 
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and also the vectors {L\ and \R) 
(L| = (1,0), \R)-- 



T- 

J- 



(8) 



(9) 



Now using ^ and ^ one can readily see that ([3]) gives 

Following [3| we define two contact operators — 
|l+)2i-i(0+|2i and bi = |1")2j-i (0~ |2i associated with 
a contact from above and below the x-axis respectively. 
Now one can calculate the probability of finding a contact 
at site 2i above or below the horizontal line as follows 



{aijn = 
{bi)n = 



{L\T'a,T'^-'\R) 
(L\T''\R) '■ 

{L\T%T'^~'\R) 
{L\T^\R) 



(10) 
(11) 



It is now easy to check that these probabilities can be 
written as 



{ai)n = 

{k)n = 



Zi{zi,0o)Zn-iizi,Z2) 

Zn(zi,Z2) 
Zi^ljzi, Z2)Zn,-i+l{00, Z2) 
Zn{zi, Z2) 



(12) 
(13) 



Using the fact that Zn{zi, Z2) — Zn{z2, zi), it can be seen 
that these probabilities have the following symmetry 



which can be understood from the symmetry of the paths. 
Let us define the average number of contacts from above 
and below the horizontal axis at position x — ^ {0 < x < 
1) as 



(15) 



Now in the thermodynamic limit n — > 00 and in different 
phases of the one-transit walk model, we find for average 
number of contacts from above cc-axis 



for zi > Z2, 
(a^) = { ^ for Z2 > zi, 

1 — X for zi = Z2 



(16) 



and for the average number of contacts from below x-axis 



1 for zi > Z2i 
for Z2 > zi, 
X for zi = Z2. 



(17) 



In the following sections we introduce a general one- 
dimensional driven-diffusive system with open bound- 
aries. We discuss the case when a single product shock 
measure can develop in the system. The steady-state 
probability distribution will be given using a matrix 
method. Finally we investigate its connections with the 
walk model. 



III. SHOCK IN A DRIVEN-DIFFUSIVE 
SYSTEM 

Consider a driven-diffusive system of classical particles 
defined on a one-dimensional lattice of length n with open 
boundaries. For simplicity we consider the case where 
each lattice site can be occupied by a single particle or it 
can be empty, although it is easy to generalize the idea to 
the systems with more than one species of particles. The 
time evolution of the probability distribution function 
of any configuration of a Markovian interacting particle 
system \P{t)) is governed by a master equation which can 
be written as a Schrodinger-like equation in imaginary 
time 



-\P{t)) = H\P{t)) 



(18) 



in which H is called the Hamiltonian [3] . The matrix el- 
ements of the Hamiltonian H are the transition rates be- 
tween different configurations. For a system with nearest- 
neighbor interactions the Hamiltonian H has the follow- 
ing general form 



7J = ^ hk,k+i + hi + hn- 



fc=i 



in which 



{ai)n{zi,Z2) = {bn-i+l)n{z2,Zi) 



(14) 



(19) 



(20) 
(21) 
(22) 
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For our two-states system T is a 2 x 2 identity matrix and 
/i is a 4 X 4 matrix describing the bulk interactions. The 
two matrices /i^'') and /i(') are both 2x2 matrices which 
describe the interactions at the boundaries. 
A shock in such a system is defined as a sharp disconti- 
nuity between a hight-density and a low-density region. 
We assume that a single product shock measure in the 
system defined as 



|fc) 



1-pi 
Pi 



1 - P2 
P2 



i — k 



(23) 



in which k {Q < k < n) is the position of the shock, can 
evolve in time generated by the Hamiltonian of the sys- 
tem. The quantities pi and p2 are the densities of the 
particles on the left-hand side and the right-hand side 
of the shock position. In [l^l the authors have studied 
the case where the shock position in has a simple 
random walk dynamics in the bulk of the system, while 
reflecting from the left and the right boundaries, by ap- 
plying some constraints on the microscopic reaction rates 
of the system. If the shock front hops to the left (right) 
with the rate 5i {5r)^ we assume that it reflects from the 
right (left) boundary with the same rate. As we will dis- 
cuss this will not alter the generality of the problem. 
The steady-state of the system it this case can be writ- 
ten as a linear superposition of the shocks ((23| . In the 
same reference the authors have also shown that the same 
steady-state can be obtained using a matrix-product ap- 
proach. According to this approach the steady-state of 
the system \P*) can be written as 



provided that one chooses the operators E and D, asso- 
ciated with the presence of a hole and a particle at each 
lattice site, as follows 

The matrix elements of the vectors \ V) and {W\, defined 
as 



{W\ ^{Wi W2) 



(26) 



besides do should only satisfy two constraints 



do—=p2~Pi , do— = (y-)(pi - P2)- (27) 

W2 Vi Or 



The normalization factor in (j24p is given by Z„ — 
{W\{D + E)"'\V). One should note that the only dif- 
ference between the representation introduced here with 
the one used in is that here we have chosen the upper 
triangular matrices instead of lower triangular matrices. 



IV. CONNECTION BETWEEN THE WALK 
MODEL AND THE DRIVEN-DIFFUSIVE 
SYSTEM 

In what follows we show that the matrix-product par- 
tition function defined in the previous section is ex- 
actly equal to the partition function of the walk model 
Z„(zi,Z2) defined in ^ provided that the fugacities of 
the contact point are chosen as 



2:1 = 1 



Z2 



Si 



(28) 



Defining C :— D + E and using a similarity transforma- 
tion we can write 

in which 

C = U-^CU , \V) = C/^V> > {W\ = {W\U. (29) 
We introduce U as 



U : 



77 



a,, 
U2 



(30) 



in which ui and U2 are assumed to be non-zero. Since 
the matrix elements of the vectors {W\ and \V) should 
only satisfy we define these elements as 



Vl 



Ul 



(24) provided that 



V2 



U2 



U2 , Wi 



1 
Ul 



, W2 



1 / 5r s 



do— = {pi - P2){ 

Ul 



) 



Now it is easy to verify that 

C^T,\V) = \R), {W\ = {L\ 



(31) 



(32) 



(33) 



which means that the partition function of the walk 
model Zn{zi = 1,Z2 = j^) is equal to the normaliza- 
tion factor of the driven-diffusive system obtained from 
the matrix-product approach in (IM)) . 
Finally we show that the density profile of the particles 
in the driven-diffusive system in terms of the matrix- 
product approach defined as 



{W\C^-'^DC''-'\V) 
{W\C'^\V) 



\<i<n 



(34) 



can be expressed in terms of the probability of finding a 
contact above or below the a;-axis defined in (fTO|) and 
pTjl . Using the same similarity transformation men- 
tioned above one finds 



{P^) 



{L\T'-'^DT"-'\R) 
(L|r«|i?) 



(35) 
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in which D = U~^DU. Now it is a straightforward step 
to show that the density profile of the particles can be 
rewritten as 

{Pi) = Pi{ai)n + P2{bi)„ (36) 

given that ([32)) is satisfied. Note that in terms of the 
fugacities zi and Z2 the phase diagram of both systems 
can be explained as follows: for zi > zi {z\ < zi) and 
in the thermodynamic limit n — > oo the walk has only 
contacts from below (above) the horizontal axis. This 
is equivalent to the fact that the shock front has a ten- 
dency to move to the left (right) which in turns means 
that the mean density of the particles in the bulk of the 
driven-diffusive system is equal to p2 (pi). One can see 
that these are confirmed by ([36]) . On the other hand on 
the coexistence line one has z\ — Z2 or in terms of the 
shock position hopping rates br — Sr- On this line the 
shock position can be anywhere on the lattice with equal 
probabilities which results in a linear density profile for 
the particles on the lattice. Using and P7)) one can 
see that gives the right density profile. 

V. CONCLUSION 

In this paper we considered two different systems: a 
one-transit walk model defined on a diagonally-rotated 



square lattice and a general driven-diffusive system in 
which a single product shock measure has a simple ran- 
dom walk dynamics. We showed that how the phase dia- 
grams of these two systems are connected by associating 
the fugacities in the walk model with the hopping rates 
of the shock front in the driven-diffusive system. In fact, 
following the idea of ^Td\ one can consider the normal- 
ization factor of a non-equilibrium system as a partition 
function of a two-dimensional equilibrium model. The 
fugacities defined in the walk model are related to the 
hopping rate of the shock front. The main difference of 
this paper with those previous is that our example be- 
longs to a larger family of driven-diffusive systems instead 
of just the ASEP as a very special example. The physical 
quantities of these two systems are directly related (for 
instance see p6p ). As a simple generalization, it can be 
shown that by redefining the fugacities of the first and 
the last steps in the walk model, one can also generate 
the partition function of the driven-diffusive systems in 
the case where the hopping rates of the shock front in 
the bulk are different from those of the boundaries. As 
another generalization to this work one can try to in- 
vestigates the driven-diffusive systems in a discrete time 
updating scheme and find the proper walk models which 
are related to them. 
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